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Abstract 

We consider systems of interacting bosons confined to one-dimensional har- 
monic traps. In the hmit of perturbatively weak two-body interactions the 
system exhibits several universal states that are exact solutions for a large 
class of two-body interactions. These states are closely related to the exact 
solutions found previously in rotating Bose-Einstein condensates. 

I. INTRODUCTION 

Bose-Einstein condensation in dilute atomic gases has received considerable attention in 
the past few years |1[ . A recent experiment reports the realization of a quasi one-dimensional 
Bose-Einstein condensate in a highly anisotropic, cigar-shaped, three-dimensional trap 
0]. Corresponding theoretical studies predict qualitatively different regimes for the con- 
densate depending on the density and the strength of the two-body interaction. For high 
densities and strong repulsive interaction, the Thomas-Fermi regime applies. Presently, this 
regime is accessed in the experiment 0. For low densities and sufficiently strong repulsive 
interactions, theory predicts a Girardeau- Tonks gas of impenetrable hard cores @-§[. A 
necessary condition for this regime is that the interaction energy per particle exceeds the 
oscillator spacing If this condition is violated, one approaches the regime of a gas with 
perturbatively weak interactions. 

In this work, we want to consider the regime of low densities and perturbatively weak 
interactions. Within this regime the interaction simply lifts the enormous degeneracy of 
the harmonically trapped A^-boson system. We will show that the Hilbert space and the 
structure of the Hamiltonian is closely related to the problem of rotational states in weakly 
interacting harmonically trapped Bose systems in two spatial dimensions . This allows 
us to transfer exact solutions from the latter problem to the present case. We will present 
universal wave functions that are exact solutions for a wide class of two-body interactions. 

This paper is divided as follows. In the following section we state the problem in first 
quantization. Next, we establish the close relationship to rotating Bose-systems. This allows 
us to present exact solutions for the one-dimensional system in analogy to the case of rotating 
systems. In the fourth section we investigate the evolution of the spectrum under change 
of a parameter and find the coexistence of exact solvability and nonintegrability. We finally 
give a summary. 
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II. HAMILTONIAN AND HILBERT SPACE 



Consider bosons in a one-dimensional harmonic trap that interact via two-body con- 
tact interactions. The Hamiltonian is 



H 



N 

■E 



+ 9 E H^i 

l<i<j<N 



(1) 



where 5^ is a couphng constant. For vanishing g = the spectrum consists of degenerate 
sets of levels that differ by multiples of the oscillator spacing. For simplicity we set the 
oscillator spacing to one and set the ground state energy to zero. Then the energies are 
simply given by integer values E = 0, 1,2, . . .. Many-body basis functions are completely 
symmetrized products (xi) . . .0„jy(xAr) of single-particle oscillator wave functions (pnix), 



0,1,2,... and fulfill E 



N 



'i=i 



E. At energy E the number of degenerate states 



n 

equals the number of partitions of E into at most integers: Hilbert space is partition 
space. Switching on the two-body interaction lifts the degeneracy between the levels. For 
sufficiently weak interaction A^l^'l -C 1 the quasi-degenerate levels are much more closely 
spaced than the oscillator spacing, and levels with different energies E do not mix. This 
is the regime we are interested in. In first order perturbation theory, the spectrum can be 
obtained by diagonalizing the interaction in the space of degenerate levels. 

We analytically treat the problem in first quantization. Let a] and dj create and annihi- 



late an excitation of boson j, i.e. d 



v^rr+T0„+i(xj), and cij 



ncpn-iixj). 



These operators fulfill the usual commutation relation [aj,d^] = 6ij. It is useful to express 
the contact interaction in terms of creation and annihilation operators: 
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We have used the Baker-Hausdorff formula in the step from the second to the third line in 
eq. (@). For perturbatively weak interactions, only the terms in the double sum with m = n 
are relevant since they conserve the total number E of excited quanta. Thus, 
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r(m + i 
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Note that a large class of two-body interactions which depend only on the two-particle 
distance can be written similarly to eq. (|^). We therefore consider generalized two-body 
interactions of the form 
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V =Y. C-mA.m, (4) 



m=0 

with 

, m 



l<k<l<N 



The coefficients depend on the details of the specific two-body interaction under consid- 
eration. One finds, e.g., that a monomial interaction of the form [xk — xiY"- yields coefficients 

(2n)!2-" 

for m < n and Cm = otherwise. 



[n — m)\ {m\y 

This result can be derived by employing the characteristic function that generates the mo- 
ments of the position variable in the harmonic oscillator. Thus, zero range interactions and 
interactions that are analytical in the two-particle distance are of the general form @ and 
(I)- 

The operator 

^ = E4«fc (6) 

k 

counts the number E of excited quanta and clearly commutes with the interaction @). In 
harmonic systems with interactions that depend only on the two-particle distances, one may 
separate the motion of the center-of-mass mode 

1 ^ 

«1 = E 4- (7) 

k=i 

The energy associated with this mode is given by 

E, = Na\ac. (8) 

This operator commutes with the interaction (^) and the energy . It has integer quantum 
numbers denoted by i?c, and = 0,1,2 . . . , E — 2, E for fixed energy E. 

III. ONE-DIMENSIONAL CASE AND YRAST LINE PROBLEM 

When stated in terms of the operators V, E and Ec presented in equations (H), (^, and 
respectively, the problem of interacting bosons confined by a one-dimensional harmonic 
potential is very similar to the Yrast line problem^ of rotating bosons in spherically sjth- 
metric two-dimensional harmonic traps [0-^. The Yrast line problem possesses analytical 
solutions for several wave functions and energies p|-[16|. We may thus transfer these solu- 
tions to the present case. Before we do so, we briefiy remind the reader of the Yrast line 
problem and its results that are relevant for us. 



^The terminology "Yrast line" is borrowed from nuclear physics. 
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The Yrast line problem consists of finding the ground states of a rotating many-body 
system as a function of total angular momentum L. Let us assume that the bosons are 
in their ground state with respect to excitations in the ^-direction. The problem thus 
becomes essentially two-dimensional. Again, one considers perturbatively weak interactions 
that simply lift the degeneracy of the harmonic trap. Single-particle wave functions are 
'iIji{z) = (27r/!)~^/^2:', where z = x + iy is the coordinate, and we have omitted the Gaussian 
in the wave function. Many-body basis states are symmetrized products of single-particle 
wave functions subject to the requirements that the single-particle angular momenta fulfill 
J2f=i = L. The Hilbert space is thus partition space and isomorphic to the Hilbert space 
of weakly interacting bosons in one- dimensional traps once we identify L = E. Let us now 
turn to operators in the Yrast line problem. Clearly, zj and ^ are creation and annihilation 
operators, respectively, when acting on the single-particle states ipi^x). Upon substitution 

a] < — > Zj, 



aj < — > 



_d_ 

dzj 



(9) 



we indeed obtain all relevant operators for the Yrast line problem |T^: Two-body interactions 
V have the form of eq. @, the total angular momentum corresponds to the operator E in 
eq. d^), and the angular momentum of the center of mass is given by the operator Ec in eq. 
(1). Thus, there is a correspondence between the Yrast line problem and interacting bosons 
in one-dimensional harmonic traps. This correspondence is based on the isomorphism of the 
Hilbert spaces and a formal identity between operators. The Yrast line wave functions 

N 

mL{zi,...,ZM)= {z.p^- Zc){zp^- z^) ...{zj,^- z^)\{e-^\^^\'' (10) 

l<Pl<P2<-<PL<N 3 = 1 

are exact solutions for a large class of two-body interactions for total angular momentum L, 
and 2 < L < TV [|13|,|l4|,|Tg . Here z^ = N-'^ J2f=i Zj denotes the center of mass. The associated 
eigenenergy is ^N{N — l)co + NL{ci + 2c2), and the coefficients cq, ci and C2 are determined 
by the specific interaction (^ under consideration. Let us transfer this important result to 
the present case. It is a peculiarity of the Yrast line problem that the creation operators, 
single-particle coordinates and single-particle wave functions are all denoted in terms of 
Zj,j = 1, . . . ,N. Clearly, the Gaussian in the wave function (p!0| ) involves only coordinates 
and no operators. Taking the correspondence between operators (P), we thus find from eq. 



(|ig) that 

N 

<^Eixi, ...,xn)= e (^ji (aL - «J) n e"^"^' (11) 

l<Pi<...<PE<Af j=l 

is an eigenfunction of interaction (^. The corresponding energy is |jl4| 



e{N, E) = ]^N{N - 1) Co + NE{ci + 2c2), (12) 
for 2 < < A^. This is the main result of this work. A direct proof of this statement may 



be obtained by repeating the derivation in ref. |jT^ and using the correspondence (P). It is 
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based on the observation that the state ([TT|) is an eigenfunction of the operators A^- The 
eigenvalues are ^N{N — 1), NE and 2NE for Aq, Ai and A2, respectively. The operators Am 
with m > 2 annihilate this state and thus have zero eigenvalue. By construction, the state 
$E is totally symmetric in the single-particle coordinates, involves E excited quanta, and 
does not excite the center-of-mass mode (^, i.e. Ec^e = 0. Note that the wave function 
(prp does not depend on the details of the interaction (^) since there is no reference to the 
coefficients Cj. It is therefore a universal solution for a large class of two-body interactions 
The corresponding eigenenergy (|l2l), however, does depend on the ffist three coefficients 
Co, ci and C2. 

Note that we may also establish a correspondence between basis wave functions of the 
Yrast line problem and the present problem. In the former system, wave functions are 
Gaussians multiplied by homogeneous polynomials of degree L that are totally symmetric 
in the single-particle coordinates Zj,j = 1,...,N. Corresponding wave functions for the 
present case are obtained by using the correspondence rule (|^) and act with the resulting 
operator on the Gaussian ground state. In this basis, both wave functions (|r^) and (|TTp 
have identical structure, i.e. identical expansion coefficients. These wave functions are also 
identical in second quantization. 

We next ask the question of whether other states in the spectrum can be written as exact 
algebraic functions. Leaving aside center-of-mass excitations, there are at most two states 
with Ec = in the spectrum when E < 6 and they are both algebraic. The representation 



of the second state is presented in ref. [14]. We believe that there are no additional universal 



states when > 6, as will be seen in the next section. 



IV. OTHER UNIVERSAL SOLUTIONS? 

Let us investigate whether there are further universal wave functions. To this purpose 
we may take an interaction (^) that depends on a parameter r and consider the spectrum 
as a function of the parameter. Avoided crossings and level repulsion would certainly in- 
dicate nonintegrability. We may further monitor the wave function structure as a function 
of the parameter and thereby search for other universal states. This requires numerical 
calculations. 

Let us consider the interaction 

W{T) = {l-r)V^a + rV^,,t (13) 

which interpolates between the contact interaction Via in one-dimensional systems and the 
interaction V^ert as r evolves from zero to one. We take V^ert = 10~'^A4 for definiteness. 
Numerical computations are most conveniently done in second quantization. Let bl^ and fe„ 
create and annihilate a boson in the single particle state \n) with energy n, respectively. 
Clearly, we have {x\n) = At fixed energy E, Hilbert space is spanned by many-body 

states \no,ni, . . . jUe) with Y,f=o n-j = N and Y.f=oj = E. The matrix elements for the 
contact interaction Via involve an integral over four oscillator functions and can be calculated 
analytically [1^. The second quantized form of the operator reads 



Vw = \Y.^^Jk^h\h]khu (14) 

ijkl 
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with 



_ 2_± /l!Tr r(t + i)r(z-t + |)r(j-t + i^ 

The operator has similar structure and matrix elements 



a. 



ijkl 



i'=max (0,m—k,l—i) 



We compute the spectrum of W{t) for < r < 1. We have fixed ii^ = 12, = 30 and 
have restricted ourselves to the states with Ec = 0. The lowest energy is set to zero. Fig.|l| 
shows several avoided crossings as the parameter r evolves, and there are no crossings. This 
indicates that the operator W{t) is nonintegrable. Inspection of the eigenstates shows that 
only the lowest energy state is independent of r. On the one hand, this finding demonstrates 
the universality of the wave function (]Tl|). On the other hand, it shows that there are no 
further states that are universal solutions. The system of bosons in one-dimensional har- 
monic traps, as well as the Yrast line problem, thus exhibits exact solvability in coexistence 
with level repulsion and nonintegrability. 

Due to their nonintegrability, the systems considered in this work differ from other exactly 
solvable many-body systems. We recall that the Calogero-Sutherland models with long- 



ranged inverse square potentials are completely integrable |]18|,|19[. Other examples are the 
one-dimensional system of zero-range interacting bosons confined to a box which is exactly 
solvable for all interaction strengths or the Tonks- Girardeau gas of hard-core bosons in 



one-dimension [^T],^. However, the present work shows that harmonically confined Bose 
systems display a few universal wave functions that are exact solutions for a large class of 
perturbatively weak two-body interactions. In this sense these fall into the class of partially 



solvable quantum many-body problems |23 



Let us finally compare the spectra of the contact interactions in one-dimensional traps 
and the Yrast line problem. To this purpose we choose V^ert = Wrast and repeat the numerical 
calculation. The contact interaction VVrast for the Yrast line problem is of identical structure 
as the operator Vid in eq. (p!4D but has matrix elements H] 



We compute the spectrum of W{t) for < r < 1. We choose = 14, = 100 and restrict 
ourselves to the states with = 0. Fig. ^ shows the result. The lowest energy is set to 
zero, and the spectra at r = and r = 1 are scaled to cover similar spectral ranges for 
display purposes. Clearly, the operators Via and VVrast yield very similar spectra and are 
in this sense close together. This suggests that one may transfer several results from the 
Yrast problem to the present case. E.g., the structure of low-lying excitations |l8|,|2^-p8| are 
expected to be similar in both problems. Though the spectrum does not undergo dramatic 
changes, we found only one universal state and a few avoided crossings as well. 
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V. SUMMARY 



We have established a close relationship between the problem of interacting bosons in 
one-dimensional harmonic traps and the Yrast state problem for interacting bosons in two- 
dimensional isotropic harmonic traps. The Hilbert spaces of both problems are isomorphic 
to each other and the important operators have identical structure. This allowed us to 
transfer several analytical results concerning eigenstatcs from the Yrast problem to the 
one-dimensional case. In particular, the one-dimensional problem exhibits universal wave 
functions that are exact solutions for a large class of two-body interactions, too. Among 
these are the lowest-lying excitations of the Bose-Einstein condensate. It is remarkable 
that one can thus learn something about rotational states from wave functions in one- 
dimensional, non-rotating systems. We observed level repulsion in a parametric Hamiltonian. 
This indicates that the system is nonintegrable. 

The author thanks George F. Bertsch for useful discussions and a careful reading of the 
manuscript. He acknowledges support as a Wigner Fellow and staff member at Oak Ridge 
National Laboratory, managed by UT-Battelle, LLC for the U.S. Department of Energy 
under contract DE-AC05-00OR22725. 
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FIG. 1. Spectrum of (1 - r) Vid + IQ-"^ t A4 as a function of r. There are several avoided 
crossings. 
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FIG. 2. Spectrum of (1 — r) Vi^ + r VVrast as a function of r. The spectrum exhibits only a few 
avoided crossings. 



10 



